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Abstract 

We discuss the first moments of the parity-even Generalized Parton Distributions (GPDs) in a 
nucleon, corresponding to six (generahzed) vector form factors. We evaluate these fundamental 
properties of baryon structure at low energies, utilizing the methods of covariant Chiral Perturbation 
Theory in the baryon sector (BChPT). Our analysis is performed at leading-one- loop order in 
BChPT, predicting both the momentum and the quark-mass dependence for the three (generalized) 
isovector and (generalized) isoscalar form factors, which are currently under investigation in lattice 
QCD analyses of baryon structure. We also study the limit of vanishing four-momentum transfer 
where the GPD-moments reduce to the well known moments of Parton Distribution Functions 
(PDFs). For the isovector moment {x)u-d our BChPT calculation predicts a new mechanism for 
chiral curvature, connecting the high values for this moment typically found in lattice QCD studies 
for large quark masses with the smaller value known from phenomenology. Likewise, we analyze the 
quark-mass dependence of the isoscalar moments in the forward limit and extract the contribution 
of quarks to the total spin of the nucleon. We close with a first glance at the momentum dependence 
of the isoscalar C-form factor of the nucleon. 



^Work supported in part by BMBF and EU-I3HP 



1 Introduction 



Understanding the structure of the nucleon arising from the underlying dynamics of quarks and 
gluons — governed by the theory of Quantum-Chromo-Dynamics (QCD) — is still one of the major 
challenges in nuclear physics. About a decade ago the concept of Generalized Parton Distributions 
(GPDs) has emerged among theorists, constituting a universal framework bringing a host of seemingly 
disparate nucleon structure observables like form factors, moments of parton distribution functions, 
etc. under one theoretical roof. For reviews of this very active field of research we refer to [1]. 

Working in twist-2 approximation, the parity-even part of the structure of the nucleon is encoded 
via two Generalized Parton Distribution functions IP{x,(,,t) and E^{x,(^,t). For a process where 
the incoming (outgoing) nucleon carries the four-momentum p^* (pi^) we define two new momentum 
variables 



The GPD variable x then corresponds to Bjorken-x, which in the parton model of the nucleon can 
be interpreted as the fraction of the total momentum of the nucleon carried by the probed quark q. 
t = denotes the total four-momentum transfer (squared) to the nucleon, whereas the "skewdness" 
variable ^ = —n ■ A/2 with n ■ p = 1 interpolates between the t- and the j;-dependence of the GPDs. 
Throughout this work we focus on the non-perturbative regime of QCD: -t << 1 GeV^, i.e. the realm 
of strongly interacting matter in our solar system making up humans, planets, etc.. We are therefore 
utilizing the methods of Chiral Effective Field Theory (ChEFT) for our analysis [2]. 

The 3-dimensional parameter space of GPDs is vast and rich in information about nucleon structure. 
The experimental program of their determination is only at the beginning at laboratories like CERN, 
Desy, JLAB, . . . [l]. However, moments of GPDs can be interpreted much easier and are connected 
to well-established hadron structure observables. E.g. the zero-th order (Mellin-) moments in the 
variable x correspond to the contribution of quark q to the well known Dirac and Pauli form factors 
Flit), F2{t) of the nucleon: 



For the case of 2 light flavors the isoscalar and isovector Dirac and Pauli form factors of the nucleon 
have been studied at low values of t at the one-loop level in Chiral Effective Field Theory, both in 
non-relativistic [3] and in covariant ^ [5] schemes. The chiral extrapolation of these form factors for 
lattice QCD data with the help of ChEFTS has been discussed in refs.[8l[9l[T0l[TT]. whereas the status 
of the experimental situation is reviewed in ref.|12j. 



^In a field theory like ChEFT the result obtained at a given order in the perturbative expansion is independent of the 
the choice of regulator applied to the UV-limit of the diagrams. (For an explicit demonstration we refer to the comparison 
between cutoff- and dimensional-regularization for the mass of the nucleon in HBChPT in ref.[6]). Nevertheless, forms 
of UV-regulators for HBChPT loop diagrams have been proposed in the literature, which effectively amount to modeling 
short-distance contributions of higher orders. The resulting chiral extrapolation functions have then been utilized to 
extend the limit of applicability of such one loop calculations towards large momentum transfer \t\ ~ 1 GeV^ and 
towards the heavy-quark limit. Here we do not advocate such an approach, but for recent work in this direction we refer 
to ref.[7]. 



A^=p^-p^; p=ip'l+pii)/2. 
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In this work we want to focus on the first moments in x of these nucleon GPDs 
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where one encounters three generalized form factors Al^^^it), i?|Q(t), C|g(t) of the nucleon for each 
quark flavor q. For the case of 2 light flavours the generalized isoscalar and isovector form factors 
have been analyzed in a series of papers at leading-one-loop order in the non-relativistic framework of 
HBChPT, starting with the pioneering analyses of Chen and Ji as well as Belitsky and Ji |13 1ll4t [T5 l ll6j. 
In this work we want to provide the first analysis of these generalized form factors utilizing the methods 
of covariant Baryon Chiral Perturbation Theory (BChPT) for 2 light flavors, pioneered in ref.[5. 
Utilizing a variant of Infrared Regularization |17j for the loop diagrams, our BChPT formalism is 
constructed in such a way that we exactly reproduc^ the corresponding HBChPT result of the same 
chiral order in the limit of small pion massecl- 

We note that (at t = 0) a, covariant BChPT calculation differs from a non-relativistic one — provided 
both are performed at the same chiral order D — by an infinite series of terms ~ {m-,r / Mq)"^ , where 
Mq denotes the mass of the nucleon (in the chiral limit), estimated to be around 890 MeV [H]. Such 
terms quickly become relevant once the pion mass niT^ takes on values larger than 140 MeV, as it 
typically occurs in present-day lattice QCD simulations of (generalized) form factors. Aside from this 
resummation property in (I/Mq)", we note that the power-counting analysis determining possible 
operators and allowed topologies for loop diagrams at a particular chiral order (see section 13. 3p is 
identical between covariant and non-relativistic frameworks. Both schemes organize a perturbative 
calculation as a power series in (p/A^)^, where p corresponds either to a small 3-momentum or to the 
mass of the pion, whereas A-,^- w 1.2 GeV denotes the scale of chiral symmetry breaking. Finally, we 
would like to note that the first moments of nucleon GPDs have also been studied in constituent quark 
models [e.g. see ref. |20]) and chiral quark soliton models {e.g. see ref.[2T]) of the nucleon, which — in 
contrast to ChEFT — can also provide dynamical insights into the short-distance structure present in 
the generalized form factors. 

This paper is organized as follows: In the next section we specify the operators with which we 
are going to obtain information on the three generalized isoscalar and the three generalized isovector 
form factors of the nucleon. In section [3] we portray the effective chiral Lagrangean required for 
the calculation, immediately followed by the sections containing our leading-one-loop results of the 
generalized form factors in the isovector (section and in the isoscalar (section [5|) matrix elements. 
A summary of the main results concludes this paper, while a few technical details regarding the 
calculation of the amplitudes in covariant BChPT are relegated to the appendices. 

2 Extracting the First Moments of GPDs 
2.1 Generalized Form Factors of the Nucleon 

In Eqs. ()4|5p of the Introduction it was shown that the first moments of nucleon GPDs are connected 
to three generalized form factors. In lattice QCD one can directly access the contribution of quark 

^We note that this property is not shared by all regularization schemes proposed for covariant BChPT during the 
past few years [18j . 

^Technically speaking, for such a comparison one has to expand the results of loop diagrams in BCHPT in terms of 
1/(IQ-k'^F^MS), where n is determined by the order of the HBChPT calculation and _FV denotes the pion decay constant. 
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flavor q to these generalized form factors of the nucleon by evaluating the matrix element [T] 



iip'klit^D ^}q\p) = u{p') 

The brackets {. . .} denote the completely symmetrized and traceless combination of all indices in an 
operator, u (u) is a Dirac spinor of the incoming (outgoing) nucleon of mass M^v, for which the quark 
matrix-element is evaluated. In ChEFT we employ the same philosophy and also extract information 
about the first moments of nucleon GPDs of Eq. ()4|5p via a calculation of the generalized form factors 
according to Eq.(l6]). 

Studying a strongly interacting system with 2 light flavors in the non-perturbative regime of QCD, 
with the methods of ChEFT one can only directly access the singlet (s) and triplet (v) contributions 
of the quarks to the the three form factors: 

iip'lQllf^ D^}q\p)u+d = u{p') 

(7) 

i{p'\Q1{fiD^yq\p)u-d = u{p') 

(8) 
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Note that the shown 2x2 unit matrix 1 and the Pauli-matrices r'*, a = 1, 2, 3 on the right hand sides 
of Eqs. (j7l8|) operate in the space of a (proton, neutron) doublet field. 

At present, not much is knowE0 yet experimentally about the momentum dependence of these 6 
form factors. The main source of information at the moment is provided by lattice QCD studies of 
these objects {e.g. see refs.[22l [231 [23]). Given that present-day lattice simulations work with quark- 
masses much larger than those realized for u and d quarks in the standard model, one also needs 
to know the quark-mass dependence of all 6 form factors, in order to extrapolate the lattice QCD 
results down to the real world of light u and d quarks. This information is also encoded in the ChEFT 
results, typically expressed in form of a pion-mass dependence of the observables under study. (The 
connection between the explicit breaking of chiral symmetry due to non-zero quark-masses and the 
resulting effective pion-mass is addressed in section [3^ ) 

The need for a chiral extrapolation of lattice QCD results for the generalized form factors of the 
nucleon leads to a further complication in the analysis: One needs to be aware that it is common 
practice in current lattice QCD analyses that the mass parameter Mn in Eqs. (|7|8p does not correspond 
to the physical mass of a nucleon, instead, it represents a (larger) nucleon mass consistent with the 
values of the quark-masses employed in the simulation. Fortunately the quark-mass dependence of 
Mn has been studied in detail in ChEFT, both in non-relativistic |6| and in covariant frameworks |19j . 
The next-to-leading one-loop chiral formulae of ref.[l9] provide a stable extrapolation functiorH up to 



''When adding on the contributions from gluons to the isoscalar quark matrix-elements of Eq.(|7]) one obtains A2'\f{t = 

0) = 1 and B2,o''(* = 0) = [M]. 

^Recent claims in the analysis of ref. |25) reporting a possible breakdown of the chiral extrapolation formula for large 
pion masses at 0{p^) level in our opinion are due to an insufficient quark-mass dependence of the vertices entering the 
HBChPT formula at that order. This issue is under discussion [18) . 
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Fn [GeV] 


Mo [GeV] 


ci [GeV-^] 


C2 [GeV-^] 


C3 [GeV-^] 


e^(lGeV) [GeV-^] 


Sm 


1.2 


0.0924 


0.889 


-0.817 


3.2 


-3.4 


1.44 






Table 1: Input values used in this work for the numerical analysis of the chiral extrapolation functions. 



effective pion-masses ~ 600 MeV. In this work we utilize the 0{p^) BChPT result (see Appendix [U]) 

MN{m^) = Mo-4cim2 + ^4^=^ _4 + + 4ci-| arccos — ^ 



128^2^2 

3ci5^m^ 



87r2F2M2 



;2M2'°HMoy'+'^^^^- 



(9) 



in order to correct for the mass effects in Eqs. (|7l8p . The coupling constants occurring in this formula 
are described in detail in ref.[K]. Possible effects of higher orders can be estimated as ©(p^) ~ 5m xf' 
where 8m could be varied within natural size estimates —3, . . . , -|-3 [18]. 

We note that the trivial, purely kinematical effect of Mjy = MjqimT^) in Eqs. (j7l8p could induce quite 
a strong quark-mass dependence into the form factors i?2 0(^5 C'2'o(t) and might even be able to mask 
any "intrinsic" quark-mass dependence in these form factors. We are reminded of the analysis of the 
Pauli form factors F2^(t) in ref.[9], where the absorption of the analogous effect into a "normalized 
magneton" even led to a different slope (!) for the isovector anomalous magnetic moment = ~ 
0) when compared to the quark mass dependence of the "not-normalized" lattice data. We therefore 
urge the readers that this effect should be taken into account in any quantitative (future) analysis of 
the quark-mass dependence of the generalized form factors B2Q{t), C2Q{t) as well. For completeness — 
based on the extensive studies of ref.[l9] — we suggest a set of couplings in tabled to be used in Eq.Q 
from which one can estimate the impact of this effect in the matrix elements of Eqs.d?! [8|). 

Finally, we note that in the forward limit t — > the generalized form factors ^20 ~ ^) 
understood as moments of the ordinary Parton Distribution Functions (PDFs) q{x), q(x) [T]: 

(x)„±rf = A'2%{t = 0)= [ dxx{q{x)+qix))^^a. (10) 

Jo 

Experimental results exist for (x) in proton- and "neutron-" targets, from which one can estimate 
the isoscalar and isovector quark contributions at the physical point [26] at a regularization scale /x. 
In this work we choose = 2 GeV for our comparisons with phenomenolog}0. In section 14.11 we 
attempt to connect the physical value for {x)u~d with recent (preliminary) lattice QCD results from 
the LHPC collaboration [21], whereas in section [5TT] we analyse the quark-mass dependence of {x)u+d 
with (quenched) lattice QCD results from the QCDSF collaboration |22j . 



^Note that this ^-dependence is not part of the ChEFT framework, as it clearly involves short-distance physics. 
However, all chiral couplings specified in section |3] carry an implicit /^-dependence (which we do not indicate), as soon 
as they are fitted to lattice QCD data or phenomenological values which do depend on this scale. 
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2.2 Generalized Form Factors of the Pion 

The first moment of a generalized parton distribution function in a pion H^[x^^^t) can be defined 
analogously to the case of a nucleon discussed above. One obtains [Tj 

l'^dxxHlix,tt) = Al{t) + {-20''Ci{t). (11) 

The two functions Al{t), C^{t) are the generalized form factors of the pion, generated by contributions 
of quark flavor q. In the forward limit one recovers the first moment of the ordinary parton distribution 
functions in a pion: 

(x)^ = Al{t = 0)= [ dxx{q{x)+q{x)). (12) 

Jo 

In the analysis of the isoscalar GPD moments of a nucleon we encounter tensor fields directly inter- 
acting with the pion cloud of the nucleon. One therefore needs to understand the relevant pion-tensor 
couplings in terms of the 2 generalized form factors ^7r(t), CT^{t). We note that the two generalized 
form factors of the pion have been analysed at one-loop level already in ref.[27] for the total sum 
of quark and gluon contributions, whereas the quark contribution to the form factors as defined in 
Eq. lllip has been the focus of the more recent work [28j . 



3 Formalism 

3.1 Leading Order Nucleon Lagrangean 

The well-known leading order Lagrangean in BChPT is given as [3] 

^^^l = \irD^ - Mo + ^7^5^^^} ^N, (13) 

with 

-D^V'Af = I^^M ~ + ^ u\d^u -^u^ {v^ + a^)u-\{v^-a^)v)^iiN-, (14) 

= iv)V^Uv}. (15) 

U = V? corresponds to a non-linear realization of the (quasi-) Goldstone boson fields, u^, d^ denote 
arbitrary isovector vector, axial- vector background fields, while v^l^'^ is the isosinglet vector background 
field. The covariant derivative V is defined as 

V^U = d^U -i {v^ + d^)U + iU {Vf, - o^) . (16) 

Finally, we note that the coupling denotes the axial-coupling of the nucleon (in the chiral limit), 
whereas Mq corresponds to the nucleon mass (in the chiral limit). 

We now want to extend this Lagrangean to the interaction between external tensor fields and a 
strongly interacting system at low energies. In this work we focus on symmetric, traceless tensor fields 
with positive parity in order to calculate the generalized form factors of the nucleon. In particular, 
we utilize the (chiral) tensor structures 

= ^ (^9tia9ul3 + 9^l|39ua - ^Siiu9ap^ X (uW^'^u ± uV^'^u^^ , 

= \ (^9^la9u|3 + 9^l|39ua - ^9fiu9aP^ v'^^ ^ • (17) 
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We note, that the chhal tensor fields transform as hV^h^ under chiral rotations, where h cor- 
responds to the standard compensator field of 2 fiavour ChPT [3], whereas V^^, is a chiral singlet. 
Chiral transformation properties for all remaining structures in our Lagrangeans are standard and 



can be found in the literature. The right- and left-handed fields V^^ 
(isovector) tensor fields of definite parity u^^, a^^, i = 1, 2, 3 via 



{R,L) 



are related to the symmetric 



T 



(18) 



while Vajs denotes the symmetric isoscalar tensor field of positive parity. 

In order to study possible interactions with external tensor fields originating from the leading order 
Lagrangean Eq. (|13p . we re- write it into the equivalent form 



,(0,1) 



(19) 



with 



(20) 



The coupling (020) been defined such that it corresponds to the chiral limit value of {x)u-d 
{{x)u+d) defined in Eq. p^ . We note that the coupling a| q is allowed to be different from unity, as we 
only sum over the u + d quark contributions in the isoscalar moments, but neglect the contributions 
from gluons. While this separation between quark- and gluon-contributions does not occur in nature, 
it can be implemented in lattice QCD analyses at a fixed renormalization scale {e.g see refs.|22 t [23 t l24j). 

While the construction of the parity-even tensor interactions with a strongly interacting system 
started from the 0{p) BChPT Lagrangean, inspection of the resulting Lagrangean Eq. (jl9p reveals 
that the leading interactions actually start out at 0{p^), as we do not assign a chiral power to any 
of the tensor fields. Furthermore, symmetries allow the addition of the parity-odd tensor interaction 
~ • We finally obtain 



(0) 



(21) 



with the coupling Aa2o corresponding to the chiral limit value of the axial quantity {Ax)u-d EH]. 

The 0{p^) part of the leading order pion-nucleon Lagrangean in the presence of external symmetric, 
traceless tensor fields with positive parity then reads 



■-tiTN 



N 



h 



^2,0 



+ 



^2,0 
4Mo 



D^ + h.c] + ...^ipN 



Dy + h.c. 



(22) 



where we have introduced V° = 5" 



(0) 



. The two new couplings H^q-, 6| q ^^'^ interpreted as the 



chiral limit values of the form factors .620(0)) -^20(0) in the limit t ^ 0. No further structures enter 
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our calculation at this ordei0. Finally, we note that the coupling 62 only allowed to exist because 
we do not sum over the gluon contributions in the isoscalar moments, otherwise this form factor is 
bound to vanish in the forward limit B2^^{t = 0) = |30j . 

3.2 Consequences for the Meson Lagrangean 

The choice of assigning the chiral power to the symmetric tensor fields V^'^''' also has the conse- 
quence that the well-known leading-order chiral Lagrangean for 2 light flavours in the meson sector 
[3T] is modifled: 



^(2) _ ^rj. 



2 



■'tTTTT 



^ . . ^V^C/t (5^- + AxlV,^-") V.C/ + x^U + xU^] . (23) 

We note that the new coupling has been defined such that it corresponds to the chiral limit value 
of {x)-,^ of Eq. p2|) . It is allowed to differ from unity because we only sum over the quark-distribution 
functions in the isoscalar channel and neglect the contributions from gluons. 

The explicit breaking of chiral symmetry via the finite quark masses is encoded via 

X = 2Bo{s + ip), (24) 

if one switches off the external pseudoscalar background field p and assigns the 2-flavour quark mass- 
matrix Mq to the scalar background field s. To the order we are working here we obtain the resulting 
pion mass m^r via 

ml = 2Bo7h + 0{Mq^), (25) 

where m = (m„ + m^) /2 and Bq corresponds to the value of the chiral condensate. The other free 
parameter at this order Fq can be identified with the value of the pion-decay constant (in the chiral 
limit). 



3.3 Power-counting with tensor fields in BChPT 

We start from the general power-counting formula of Baryon ChPT: 

D = 2L + l + Y,{d-2)N!' + Y,id-l)Nf''. (26) 
d d 

D denotes the chiral dimension p^ of a particular Feynman Diagram, L counts the number of loops in 
the diagram, whereas the variables N^'^' count the number of vertices of chiral dimension d from 
the pion (M) and pion-Nucleon (MB) Lagrangeans. 

To leading order we only have the tree level contributions from the p'^-Lagrangean of Eq. (|21|) with 
L = 0, = 0, N^'^^ = 1, resulting in D = 0. 

At next-to-leading order D = 1 we find additional tree level contributions from the p^ Lagrangean 
of Eq.(l22]) with L = 0, N^^ = 0, N^^^ = 1. 

The first loop contributions enter at D = 2 with L = 1, = 0, A'^^^^ = 1 plus possible con- 
tributions from N^^. The corresponding diagrams are shown in Figdl Diagram (e) in that Figure 
represents (loop) corrections from the nucleon Z-factor (given in appendix [C]) which at this order 
renormalize only the tree-level tensor couplings of the p'^-Lagrangean. Note that there is an additional 
possibility of obtaining D = 2 contributions via L = 0, ''^ = 0, N2^^ = 1, corresponding to further 
tree level contributions discussed in the next subsection. 

^We only show those terms where the tensor fields couple at tree level without simultaneous emission of pious, photons, 
etc., as these are the relevant terms for our 0{p^) calculation of the form factors according to the power-counting analysis 
of subsection 13.31 
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d) 



Z^ 



e) 



Figure 1: Loop diagrams contributing to the first moments of the GPDs of a nucleon at leading-one- 
loop order in BChPT. The solid and dashed lines represent nucleon and pion propagators respectively. 
The solid dot denotes a coupling to a tensor field from the 0{p^) Lagrangean of Eq. (j2ip . Aside from the 
trivial unity-contribution, wavefunction renormalization to the couplings of the order Lagrangean 
only start to contribute at next-to-leading one loop order. 




Figure 2: (Isoscalar) tensor field coupling to the pion cloud of the nucleon. This process only starts 
to contribute at next-to-leading one-loop order in BChPT. 



In this first paper we will stop with our analysis of the generalized form factors at the D = 2, i.e. 
0{p'^) level, corresponding to a (covariant) leading-one-loop calculation. The next-to-leading one-loop 
effects of D = 3 will be postponed to a later communication [32] . 

The (perhaps) surprising finding of this powercounting analysis is the observation that the tensor 
coupling to the pion field controlled by the coupling in Eq. (j23p does not contribute at leading- 
one-loop order! It only starts to enter at L> = 3 via L = 1, A^2^ = 1 ^^'^ Nf^^ = 1 or 2. (The 
corresponding diagram for N^^^ = 2 is shown in FigJJl while the not-shown diagram for N^^^ = 1 is 
expected to sum to zero due to isospin-symmetry.) We therefore note that the generalized form factors 
of the nucleon behave quite different from the standard Dirac and Pauli form factors of the nucleon, 
where the pion-cloud interactions with the external source a la Figl2] are part of the leading-one-loop 
result and play a prominent role in the final result. We will discuss the impact of this particular D = 3 
contribution further in section [531 when we try to estimate the possible size of higher order corrections 
to our 0{p^) analysis. 

We now move on to a discussion of the tensor interactions in the 0{p^) Lagrangean which are 
present at D = 2 according to our analysis. 
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3.4 Next-to-leading Order Nucleon Lagrangean 



At next-to-leading order the covariant BChPT Lagrangean for 2-flavour QCD contains 7 indepen- 
dent terms in the presence of general scalar, pseudoscalar, vector and axial-vector background fields, 
governed by the coupling^ ci , . . . C7 |33j . Extending this scenario to symmetric and traceless tensor 
background fields with positive parity, the symmetries allow the construction of six additional terms 
to describe the interaction at tree-level to this ordeJfl: 



4% =^i^N[Tr{x+)V+l^iD'^ + /i.cjvTV 

+ ^^^{^''^^' + ^-"j^^ ^'■'^^ 
+ ... 



with x+ = u^x^^ + y-X^"^- The physics behind these couplings Cj, i = 8, ... 13 with respect to the 
generalized form factors of the nucleon is quite simple: cg (cg) governs the leading quark-mass insertion 
in {x)u-d {{x)u+d), whereas the couplings cio, cn give the values of the generalized form factors 
^2,0(0)) ^2^(0) in the double limit t ^ 0, 171-,^ — > 0. We can therefore denote them as 02^, c^ o- 
Finally, the couplings C12 , C13 parametrize the contributions of short- distance physics to the slopes of 
the generalized form factors ^20(^)5 ^2 o(*) (™ chiral limit). Note that the operator controlled by 
the coupling cg is not allowed to exist when we add the gluon-contributions on the left hand side of 
Eq.dZD [M]. 

After laying down the necessary effective Lagrangeans for our calculation, we are now proceeding to 
the results of our calculation. 



*The couplings eg and cr axe often denoted as k° and in the literature, as they constitute the leading terms in the 
isovector and isoscalar Pauli form factor of the nucleon [4]. 

^More couplings will appear when one wishes to extend this to more general scenarios. 
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4 Generalized isovector Form Factors in 0{p^) BChPT 



4.1 Moments of the isovector GPDs at t = 

In this subsection we present our results for the generahzed isovector form factors of the nucleon at 
t = 0. For the PDF-moment A^^it = 0) we obtain to ©(p^) in BChPT 



—Ja^, log -T71 + 9\^=^= 14 - 8-^ + -I arccos 
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+4m^^ + 0(p='). (28) 

Many of the parameters in this expression are weh known from analyses of chiral extrapolation func- 
tions. Numerical estimates for their chiral limit values can be found in table [H Furthermore, in a first 
fit to lattice data we constrain the coupling Aog q from the phenomenological value of {Ax)^^' w 0.21 
via [29] 

(Ax)„_d = Aalo + 0{ml) (29) 

and perform a 2-parameter fit with the couplings O20) Cg'^^(lGeV) at the regularization scale A = 1 
GeV. We fit to the (preliminary) LHPC data for this quantity as given in ref.[24], including lattice data 
up to effective pion masses of ~ 600 MeV. The resulting values for the fit parameters together 
with their statistical errors are given in table [2] and the resulting chiral extrapolation function is 
shown as the solid line in FigjSl We note that the extrapolation curve tends towards smaller values 
for small quark-masses, but does not quite reach the phenomenological value at the physical point, 
which was not included in the fit. We have therefore tested what might happen if we estimate possible 
corrections to the solid curve arising from higher orders. From dimensional analysis we know that the 

AjMo 

Constraining 6a between value^^ — 1, . . . , +1 and repeating the fit with this uncertainty term included 
leads to the grey band indicated in FigjSl Reassuringly, the phenomenological value for {x)u-d lies 
well within that band of possible next-order corrections, giving us no indication that something may 
be inconsistent with the large values for {x)u-d typically found in lattice QCD simulations for large 
quark-masses. The resulting values for the couplings Ogg; Cg^^(lGeV) of Fit I are also well within 
expectations. 

We note that the mechanism of the downward-bending at small quark- masses in A2 g (0) found 
in Eq. (j28p is quite different from what has been discussed so far in the literature within the non- 
relativistic HBChPT framework {e.g. see ref.j35]). In order to demonstrate this we truncate Eq.(|28p 



leading chiral contribution to {x)u-d beyond our calculation takes the form 0{p^) ~ 6a J^It + •••• 



The natural size of all couplings in the observables considered here is below 1, as all coupling estimates in this section 
refer to a moment of a parton distribution itself normalized to unity. For this estimate we assume that the observable 
under consideration has a well behaved chiral expansion, with the dominant contribution provided by the leading term 
in the chiral expansion. This expectation is confirmed by the fit values of tables [2] and [3] extracted for the investigated 
coupling constants. 
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Figure 3: "Fit I" of the 0{p^) result of Eq.l^ to the (preUminary) LHPC lattice data of ref.[2l]. 
The corresponding parameters are given in Table [2j Note that the phenomenological value at physical 
pion mass was not included in the fit. The grey band shown indicates the estimate of possible 0{p^) 
corrections as discussed in the text. 





Fit I (4 points - 2 parameter) 


Fit II (6+1 points - 3 parameter) 


"2,0 

cg(lGeV) 


0.157 ± 0.006 
0.210 (fixed) 
-0.283 ± 0.011 


0.141 ± 0.0057 
0.144 ± 0.034 
-0.213 ± 0.03 



Table 2: Values of the couplings resulting from the two fits to the LHPC lattice data for {x)u-d- 
The errors shown are only statistical and do neither include uncertainties from possible higher order 
corrections in ChEFT nor from systematic uncertainties connected with the lattice simulation. 



at l/iWir'^F^Mo) to obtain the exact ©(p^) HBChPT result of refs.[36l [37]: 



^lom^HBChPT = «2,o{ 1 -(^^(25U (35^ + 1) log 9')}+4m2 



Mi 



As already stated in the Introduction, the covariant BChPT scheme used in this work is able to repro- 
duce exactly the corresponding non-relativistic HBChPT result at the same order by the appropriate 
truncation in l/(167r^F^Mo). 

Fit I is certainly constricted by the assumption that we use the physical value of (Ax)^^' ~ 0.21 
for the coupling Aa2 05 which presumably takes a value in the chiral limit a bit smaller than the 
phenomenological value at the physical point [29]. Furthermore, in order to also compare the O(p^) 
HBChPT result of Eq.^ with the 0{p^) covariant BChPT result of Eq.^ we perfor m a second 
fit: We fit the covariant expression for {x)u~d of Eq. ()28p again to the LHPC lattice data, this time, 
however, we constrain the coupling Aa'2 q in such a way, that the resulting chiral extrapolation curve 

reproduces the phenomenological value of {x)^^' = 0.160 ± 0.006 [26] exactly for physical quark 
masses. The parameter values for this Fit II are again given in table [21 whereas the resulting chiral 
extrapolation curve of the covariant 0{p'^) expression of Eq. (j28p is shown as the solid line in Fig HI 



11 



First, we would like to emphasize that the curve looks very reasonable, connecting the physical point 
with the lattice data of the LHPC collaboration in a smooth fashion. Second, we note that the 
resulting values for the couplings 02 qi ^cl2,o underlying this curve are very reassuring, indicating that 
both {x)u-d and {Ax)u-d are slightly smaller in the chiral limit than at the physical point! Likewise, 

(r) 

the unknown quark-mass insertion Cg contributes in a strength just as expected from natural sizes 
estimates. For the comparison with HBChPT we now utilize the very same value£3 for 02 q and c^^ 
of Fit II as given in table [2j The resulting curve based on the 0{p'^) HBChPT formula of Eq. pop 
is shown as the dashed curve in FigUJ One observes that this leading-one-loop HBChPT expression 
agrees with the covariant result between the chiral limit and the physical point, but is not able to 
extrapolate on towards the lattice dat J^. However, we note that all our numerical comparisons with 
HBChPT results shown in the figures of section H] and [5] are based on the assumption that the D=2 
fit values found in tables [2] and [3] are already reliable estimates of the true, correct values of these 
couplings in low energy QCD. Clearly, the shown HBChPT curves might have to be revised if future 
D=3 analyses [32] lead to substantially different numerical values for these couplings. The true range 
of applicability of HBChPT versus covariant BChPT can only be determined, once the stability of 
the employed couplings is guaranteed. A study of higher order effects is therefore essential also in 
this respect. Ideally we would like to reanalyse the results of Fit II by first fixing the couplings from 
a fit of the HBChPT formula of Eq. (l30p and then study the resulting 0{p'^) BChPT curve for this 
observable. However, a fit to the present set of lattice data shown in FigJS] leads to a curve that is 
not compatible with the 2 lightest lattice points and lies significantly below the grey band shown in 
Fig 13] — even down to the chiral limit. We will study this issue further in ref. [32| . 

Keeping these caveats in mind, at this point we conclude that the smooth extrapolation behaviour 
of the covariant 0{p'^) BChPT expression for {x)u-d of Eq. ()28p between the chiral limit and the region 

of present lattice QCD data is due to an infinite tower of (^jif^^ terms. According to our analysis 

the chiral curvature resulting from the logarithm of Eq. (|30p governing the leading-non-analytic quark- 
mass behaviour of this moment is not responsible for the rising behaviour of the chiral extrapolation 
function, as had been hypothesized in ref. [35] • It will be interesting to see how well our 0{p'^) BChPT 
extrapolation formula for {x)u-d will perform, once the new lattice QCD data from dynamical simu- 
lations at small quark masses for moments of GPDs become available |39| . 

Finally, we are discussing the 0{p'^) BChPT results at t = for the remaining two generalized isovector 
form factors of the nucleon at twist-2 level. One obtains 



^ log ^ , 5 - 5^ + arccos \ + O(p^), (31) 



The covariant BChPT approach appUed in this work guarantees that the chiral limit value for an observables is 
mandated to be the same in both frameworks. For this specific case this implies that there is one unique set of numerical 
values to be determined for the couplings 02.0 and c>^\ regardless of the ChEFT scheme employed. 

We note that this behaviour of the 2-flavour HBChPT result is completely analogous to the corresponding leading- 
one-loop HBChPT expressions for the axial coupling constant of the nucleon [38] and also for the anomalous magnetic 
moments of the nucleon [8]. It appears to be a pattern that such leading HBChPT extrapolation formulae only describe 
the quark-mass dependence between the chiral limit and the physical point. If one wants to push the chiral extrapolation 
to larger quark-masses, more m^-dependencies seem to be required, which from the viewpoint of the HBChPT framework 
are "higher order". 
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Figure 4: "Fit 11" of the 0{p^) BChPT result of Eq.([28]) to the (preHminary) LHPC lattice data of 
ref.[2^ and to the physical point (solid line). The corresponding fit-parameters are given in table [2l 
The dashed curve shown corresponds to 0{p'^) result in the HBChPT truncation (see Eq. (j30p ). The 
shaded area indicates the region where one does not expect that ChEFT can provide a trustworthy 
chiral extrapolation function due to the large pion-masses involved, albeit the covariant 0{p'^) result 
does not signal any breakdown in this region. 
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ml 
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Mi 



+ 



Ml 



arccos 



m 



2Mo) 



+ 0(p3). (32) 



As one can easily observe, one encounters plenty of non-analytic terms and even chiral logarithms 
in these covariant Oip^) BChPT results. However, we note that in the HBChPT limit to the same 
chiral order 0{p'^) one would only obtain the chiral limit values 62 qj C2Q-> nothing else [H]. E.g. the 
chiral logarithms calculated in ref.[I5 fo^' -^2 o(* — 0)' ^2 o(* — 0) would only show up in a full 0{p^), 
respectively 0{p'^) BChPT calculation of the generalized form factors. From the viewpoint of power- 
counting in BChPT they are to be considered part of the higher order corrections to the full OijP') 
results given in Eqs. (131|32p . 

Unfortunately, at this point no information from experiment exists for these two structure quantities 
of the nucleon. From phenomenology one would expect that B'!^Q[t = 0) has a "large" positive value 
at rriTT = 140 MeV, as it corresponds to the next-higher moment of the isovector Pauli form factor 
F2 (i = 0) = = 3.7 n.m. (Compare Eq.(l5|) and Eq.(l3]) at ^ = 0). Lattice QCD analyses seem to 
support this expectation [221 [23l [21! • 1^ contrast, the value of C2o(* = 0) cannot be estimated from 
information known about nucleon form factors. State-of-the-art lattice QCD analyses {e.g. see ref.|24]) 
suggest that it is consistent with zero. In FigOwe have indicated how the corresponding extrapolation 
curves based upon this information might look like. We are looking forward to a combined analysis 
of Eqs (j31|32p and the new lattice QCD data from dynamical simulations at small quark-masses [39] 
extrapolated to t = 0, which will shed new light onto this new domain of nucleon structure physics! 
As a caveat we note that in both form factors we would see a dominant influence of the quark-mass 
dependence stemming from the kinematical factor Mpf{mTr) in Eq.dS]), z/ their corresponding chiral 
limit values &20' "^20 7^ ^- ^ further observation is that the uncertainties connected with possible 
higher order corrections from 0{p'^) could already become substantial for pion masses around 300 
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Figure 5: Quark-mass dependence of the isovector moments -Bg = 0) and C|Q(f = 0). In -B|o(* ~ ^) 
we have varied the (unknown) chiral hmit value 62 between and +0.5, as lattice analyses [221 [231 [23] 
suggest that this moment has a large positive value. For the chiral limit value of C2o{t = 0) we 
have chosen the value 0, as preliminary lattice QCD analyses suggest that this moment is consistent 
with zero [24]. The grey bands shown indicate the size of possible higher order corrections to these 
0{p'^) results. 



MeV. For both quantities they can be estimated vi4__l 0{p ) ~ 5b,c 



m^Afiv(m^) 



where — 1 < 



,c 



< 1. 



In order to ultimately test the stability of the results in Eqs. ()31|32p it will be very useful to extend 
this analysis to next-to-leading one- loop order [32] . 



4.2 The slopes of the generahzed isovector Form Factors 

In order to discuss the generalized isovector form factors ^2o(*)' -^2o(0' ^2o(^) non-zero values of 
t, we first want to analyze their slopes px, defined via 



X = A,B,C. 



(33) 



To 0(p2) in BChPT we find 



Pi 



C12 _ 0'2,o9a ml 

6(47rF^)2 {AM^ - ml) 



mi / 39 , mi 



26 + 8 log 



mt 



90 



fi 2 
-3^1og^ 

2 4 

mt, ml. 







mi 
Ml 



M^ 



30 + 32 log 



arccos 



2Mn 



+ 



m„ 



AlMo 



Ai 



(34) 



^■^Note, that the momentum in Eas. (|7l8p counts itself as an object of order p leading to the observation that a 
D = 3 evaluation of the currents Eas. (|7l8| ) leads to a leading quark-mass dependence of m^~^ in i32''o'''(0). 
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) 



(36) 



(35) 



) 



The parameters (5^, S^^, 6^ are not part of the covariant 0{p'^) result. They are only given to indicate 
the size of possible higher order corrections from 0{p^). A numerical analysis of the formulae given 
above suggests that the size of pion-cloud contributions to the slopes of the generalized isovector 
form factors is very small! The physics governing the size of these objects seems to be hidden in 
the counterterm contributions Ci2, S^^, 6^ which dominate numerically when assuming natural size 
estimates — 1 < ci2, S^^, 6^ < +1. We note that this situation reminds us of the isoscalar Dirac and 
Pauli form factors of the nucleon -Ff (t), -^1 (^)' where the t-dependence in SU(2) ChEFT calculations 
is also dominated by counterterms {e.g. see the discussion in ref.[3]). 



Finally, truncating the covariant results of Eqs . (13411 36]) at l/(167r^F^Mo) we obtain the (trivial) 
HBChPT expressions to 0{p'^): 



The non-zero slope results found in the HBChPT calculations of refs.[15] [T6] are therefore of higher 
order from the point of view of our powercounting. Most of them can already be added systematically 
to our covariant ©(p^) results of Eas. pill36]) aX. 0{p^) [32]. 



4.3 Generalized isovector Form Factors of the Nucleon 

In this subsection we present the full t-dependence of the (generalized) isovector form factors of a 
nucleon to 0{p'^) in BChPT. We note that for all three form factors at this order only diagram c) 
of Figd] gives a non-zero contribution at finite values of t — the resulting expressions at this order are 
therefore quite simple: 




+ (!?(/), 
+ O(/). 




(37) 



(38) 
(39) 



A. 



it) 



A 



(0) + 



Floit) + ^t + 0{p^) 



(40) 



1927r2F2 
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with 



t ; — =: arccos 



MoV4M^ - ml \2Mo, 

lOM^ + (I7?n2 + QOM^) Af^ - 4ml {^1 + 15M^) 



\ J 2mlMl 



I Vm^xUm'^ - ml 



+ 



^Vlm%Ml\ arccos ( | + 



2 ( Mq^ - 



MlM^ 



+ 



2Mo2 - 



yim%M^ + 2m2 (4m2 - QM^) 



2M* + Sm^ (3m2 + 4Mo2) 



-4m^ [ml + OMq^) + UmlM^ 



log I IT- I + 



2mi 



Mo ) M^M^ 



mt (-12M^ + AM'^M^ 



log 



Mo 



(41) 

and = Mq + (u^ - i) L Note that F2^Q{t = 0) = 0, whereas ^2,0 (0) heen discussed in section 

14. 1[ A conservative estimate for the size of possible higher order corrections indicated in Eq. (j40 '' 

„. ^3 



can 



3 

be obtained via 0{p^) ~ 6a + ^\ a^Mo ^' ^ discussed in the two previous subsections. For the 
remaining two (generalized) isovector form factors we obtain 



"2,0 



Mo 



+ 



-GmlM^M"^ + em^Mo^ (m^ - Sm^ M^ ^ ^ 



6ml 



4M2 



m 



.2 L 



mt - 5mlM^ + 5M^ 



arccos 



V2M 



AlMo 



+ Ob 



(42) 



cut) 



^2,0 ~ r 



Mo 



487r2F2 



i MS 



+6m^Mo2M2 - em^Mo^ [ml - 2M^j log 



m^ 



+ 



+5c 



dinlM^ 



4M2 



mt - AmlM"^ + 2M^ 



arccos 



A2Mo 



^ MAr(m^ 



A2Mo 



2 ( M|5^ - M^ ) M^ - Sm^Mn^M 



V2M/ 



(43) 



The parameters 6b, 6c, (5^, 6^ have been inserted to study the possible corrections of higher orders 
to our 0{p'^) results. Varying these parameters between -1 and +1 we conclude that a full 0{p'^) 
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calculation is required, before we want to make any strong claims regarding the t-dependence of -B|Q(t) 
and C2 o(i) beyond the linear t-dependence discussed in the previous subsection. On the other hand, 
the predicted t-dependence of the form factor Al^oit) of Eq. (j40p appears to be more reliable at this 
order, as possible higher order contributions only affect terms beyond the leading linear dependence 



in t. 



5 Generalized isoscalar Form Factors in 0{p^) BChPT 

5.1 Moments of the isoscalar GPDs at t = 

To 0{p'^) in 2-flavor covariant BChPT the only non-zero loop contributions to the isoscalar moment 
^2 o(i = 0) (see Eq.([TO|)) arise from diagrams c) and e) in Fig{TJ One obtains 



^lo(O) 



{x)u+d 



+ 



mt 



2 2 



rr4 r< 



arccos | :r^^ I + 0{p^). 



V2Moy 



(44) 



Eq. (j44p should provide a similarly successful chiral extrapolation function for {x)u+d as the covariant 
0{p^) BChPT result ofEq.^ did for the LHPC lattice data for {x)u-d in sectionSTl The uncertainty 

3 

arising from higher orders can be estimated to scale as 0{p^) ~ 5^ a^Mq ' where 5^ should be a number 
between — 1,...,-|-1, according to natural size estimates. We note that the coupling Aa^'oi which 
played an essential role in the chiral extrapolation function of {x)^^^, is not present in the quark-mass 
dependence of the isoscalar moment {x)u+d- The resulting chiral extrapolation function is therefore 
presumably quite different from the one in the isovector channel. The absence of a chiral logarithm 
~ logm^ in {x)u+d (compare Eq. (f30|) and Eq. (jl^ ) presumably will only lead to a difference in the 
chiral extrapolation functions between the isovector and the isoscalar moment for values m.,^ < 140 
MeV. 

Note that from Eq.dM]) in the limit I / {IGTr"^ Mq) ^ we reproduce the leading HBChPT result 
for {x)u^(i of ref.[36], which found a complete cancellation of the non-analytic quark-mass dependent 
terms in this channel (in the static limit): 



2 eg 



^hmUchPT = ai,o + 4m^^ + 0(l/(167r^F>/o)). 



(45) 



The coupling cg is therefore scale independent (in dimensional regularization) and constitutes the 
leading correction to the chiral limit value a| g of {x)u+d- 

At t = we also find non-trivial results for the two other generalized isoscalar form factors of the 
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nucleon to this order in the calculation 

+4 log 4 - . - 54 + 41 -ccos ^ 



..1^=^, (46) 



Note that additional chiral logarithms appear in the HBChPT calculations of refs. |15^ I16| . e.g. a 
term h2Qml\ogm.,^ entering the analogue of Eq. (j46p . According to the chiral power-counting formula 
Eq. (j26p such a term only starts to appear at ~ 0{p^), and is therefore not present in our 0{p'^) 
calculation. We have added the parameters (5^, 5^ "by hand" to our covariant 0{p'^) results Eqs. (j46p 
and (|47p in order to indicate possible effects of higher order {i.e. 0{p^)) corrections. Presumably they 
take on values — 1, . . . , +1. Note that our 0{p^) BChPT prediction for o(0) is strongly affected by 
possible corrections from higher orders. This is due to the fact that one only receives a non-zero result 
for this form factor starting at 0{p^) — the order we are working in. Eq. (j47p should therefore only be 
considered to provide a rough estimate for the quark mass dependence of this form factor at i = 0. 
For a true quantitative analysis of its chiral extrapolation behaviour the complet^^ 0{p^) corrections 
should first be added [32] . 

Unfortunately, the (published) lattice QCD data base for {x)u+d and -BfoC^ ~^ 0) is quite sparse for 
small values of the pion mass. In order to obtain a rough estimate of the chiral extrapolation functions 
resulting from Eqs. (|44l46p we have to resort to (quenchedj^ data of the QCDSF collaboration in 
ref.[22]. Performing a combined fit of Eqs. (l44|46p to the lattice data shown in Fig l6] and including 
the phenomenological value of {x)u+d ~ 0.54 [26] we obtain the two solid curves shown in Figl6l The 
resulting parameters of the fit are given in table [3l Interestingly, despite the large quark-masses and 
the huge error bars in the data of ref.[22] we obtain reasonable chiral extrapolation curves with natural 
size couplings. The analysis of the QCDSF data in combination with the physical value for {x)u+d 
suggests that the chiral limit value of this isoscalar PDF-moment is smaller than the value at the 
physical point, leading to a monotonically rising chiral extrapolation function as shown in the left 
panel of FigEl It will be very interesting to observe what values an upcoming analysis of the new 
LHPC data will find [39] for this moment, as a glance at the preliminary results shown in ref.[23] seems 
to indicate that for effective pion masses around 300 MeV LHPC finds values for {x)u+d which are 
lower than the physical point. If this observation can be confirmed in the final analysis of the data [39| 
one would have to conclude that the chiral extrapolation behaviour of {x)u+d is truly extraordinary, 
displaying a minimum for effective pion masses near 300-400 MeV before rising again to larger values 
for large quark masses. 

^''The most prominent correction from O(p^) arises from the triangle diagram of Fig[2]and is given in appendix ID. 21 
via AC^.o.(^ ~ 0, m^). We note, however, that this is not the only correction at the next order. 

^^The quark-masses employed in ref. [22j are so large, that one does not expect to find differences between quenched 
and dynamical simulations. See e.g. the discussion ref.[8]. Note that we are utilizing the lattice data of ref. |22) with the 
scale set by ro, as we consider the alternative way of scale-setting (via a /ineor extrapolation to the physical mass of the 
nucleon) also discussed in ref. [22] to be obsolete in the light of the detailed chiral extrapolation studies of ref. |19). We 
further note that in the simulation results of ref. [22) all contributions from disconnected-diagrams have been neglected, 
which gives rise to an additional unknown systematic uncertainty in the lattice data. 



18 



0.65 



0.6 



A 0.55 
V 



0.5 



0.45, 



- ' 1 ' ; ' 

/ 

/ 

/ 


1 1 1 - 


/ 

/ 






/ 








/ ■ 

/ ^ 








/ 












1,1,1 


1 



0.2 



0.4 0.6 
m [GeV] 



0.8 




0.4 0.6 
m [GeV] 



Figure 6: Solid lines: The quark mass dependence of the 0{p'^) BChPT result for {x)u+d of Eq. (j44p 
and B2 (0) of Eq. (j46|) . The three unknown parameters resulting from a combined fit to the shown 
QCDSF data of ref.[22] and to the phenomenological value for {x)u+d ES] are given in table [3l At 
the physical point we obtain .B| q = —0.056 it 0.016, with the shown error only being statistical. Due 
to the poor data situation with its unknown systematic errors, however, these results should only be 
considered as a rough estimate of the true quark-mass dependence. The dashed line corresponds to 
the respective HBChPT result at this order. 



"2,0 


°2.0 


eg 




«-"-(^ = 2GeV) 


0.527 ±0.007 


-0.103 ±0.016 


0.147 ±0.002 


(fixed) 


0.538 ± 0.012 (fixed) 



Table 3: The values for the three tensor coupling constants entering ^2,0 (0 ^'^d -^2,0 (^) order 
as extracted from a combined fit to the lattice points for A|q(0) and -BloC'-') 122] shown in Figl6] 
and to the physical point of ^|o(* ~ 0,171.,^ = 0.14 GeV) = (x)^^" . Note that we have obtained a 
small negative value for 62 ■ '^^^ indicated errors are only statistical and do not reflect the (much 
larger) systematic uncertainties connected with a fit to the (quenched) lattice data of ref.[25] which 
e.g. neglect all contributions from disconnected diagrams. 



As a second observation, we would like to note that the value for the generalized form factor 
i?| Q (t = 0) could take on a small negative value at the physical point according to the right panel of 
FiglHl albeit with a large uncertainty due to the poor data situation. Because of the small (negative!) 
value of the isoscalar Pauli form factor i^l = 0) = = —0.12 n.m., it is somewhat expected that the 
next-higher momenlF^ will yield a value close to zero. However, Figl6] now opens the possibility that 
^2o(^ = 0) ~ —0.06 might be as large as 50% of its i^lC^ ~ ^) analogue. It will be very interesting 
to observe whether this feature can be reproduced when the new data of QCDSF and LHPC [39] at 
small pion mass are analyzed with the help of our 0{p^) BChPT formulae Eqs. (|44l46p . 
Finally, we note again that the true range of applicability of HBChPT versus covariant BChPT (see 
e.g. the left panel of Fig.®) can only be determined, once the stability of the employed couplings 
is guaranteed, see the similar discussion for {x)u-d in section 14.11 A study of higher order effects is 
therefore essential also in this respect. 



"We can set f = in Eq.(|3]) and Eq.((5ll for such a comparison. 
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5.2 The contribution of u— and d— quarks to the spin of the nucleon 

In the past few years a lot of interest in generalized isoscalar form factors of the nucleon has focused 
on the values of j4| q, -B| q at the point t = 0, as one can determine the contribution of quarks to the 
total spin of the nucleon via these two structures 



To 0{p^) in BChPT we find 



Ju+d = J[^i,o(i = 0)+i?|,o(i = 0)]. (48) 



'Ju+d — o 1 '^2,0 + ''2,0 77 r 



2 r^'*^ ' "^'^ Mo ' (4^F^)2 
2 / ^2 



6— -77 + — — r arccos ■ 



V4M^ -ml\ M2 MoV 2Mo 

Note that despite the plethora of non-analytic quark-mass dependent terms contained in the 0{p'^) 
BChPT result of Eq. ()49p . the two chiral logarithms calculated in ref.[13] within the HBChPT frame- 
work are not yet contained in our result. Both terms (~ a(g)7r, ^(ij)7v ™ the notation of ref.[13]) are 
part of the complete 0{p^) result according to our power-counting and will appear in the calculation 
of the next ordei0. We further note that the two logarithms of ref.[13] are UV-divergent and are 
accompanied by a counterterm, whereas the O(p^) BChPT result of Eq. (j49p happens to be UV-finite 
to the order we are working in. In ref.[T3] the authors also reported that the two chiral logarithms (of 
0{p^)) which they describe presumably are canceled numerically by pion-cloud contributions around 
an intermediate Delta(1232) state. We can confirm that this possibility exists, as the described Delta 
contributions also start at 0{p^), assuming a power-counting where the nucleon-Delta mass difference 
is counted as parameter ~ p^ (in the chiral limit) (see ref.@0| for details). 

Utilizing the 0{p'^) BChPT result of Eq. (j49p and the fit parameters of table [3] we obtain a a first 
estimate for the contribution of u— and d— quarks to the spin of a nucleon: 

Ju+d{t = 0, = 0.14 GeV) w 0.24 ± 0.05, (50) 

which is only 50% of the total spin of the nucleon! Note that we can only give a rough estimate for the 
error in this determination via assigning a typical size of the input lattice error bars shown in Figl6] 
to the extrapolated result, as we are assuming that the true error is dominated by systematic errors 
in the lattice input to our analysis. (The size of the statistical error read off from the fit of table [3] is 
±0.01 and therefore negligible. The small value of this statistical error is of course heavily influenced 
by the error assigned to the phenomenological value of {x)^^^ given in table [3]) 

Based on the same input (uncertainties) we can also predict the quark-mass dependence of Ju+d- 
The result is displayed as the solid line in FiglTl Note that in contrast to the analysis given in ref . |22j . 
we do not obtain a flat chiral extrapolation function between the lattice data and the physical point. 
The Oip^) BChPT analysis suggests that the value at the physical point lies lower than the values 
obtained in the QCDSF simulation at large quark-masses. However, given the large values of quark- 
masses and the sizable error-bars in the (quenched) simulation of ref. [22], our extrapolated value for 
Ju+d obviously can only give a first estimate of the true result. We are therefore looking forward to 
the application of our formula Eq. ()49p to the new (fully dynamical) lattice QCD results by QCDSF 
and LHPC obtained at lower quark-masses with improved statistics [39j . Furthermore, the influence of 
possible corrections from 0{p^) in Eq. (j49p has to be analyzed, before one can obtain a high precision 
determination of the quark-contribution to the spin of the nucleon from a combined effort of lattice 
QCD and ChEFT [32]. 



''The contribution ~ ffl(q),r is already contained in the function ABh.oit ~ 0,m^) discussed in subsection 15.31 
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Figure 7: The contribution of u + d quarks to the spin of a nucleon as function of the effective pion 
mass. The C(p^) BChPT result shown as the solid line is a prediction of Eq. (j49p which utilizes the 
fit-parameters of table [3l For comparison we have also plotted simulation data from QCDSF [22] in 
the figure. At the physical point one can read off Ju+d ~ 0.24. (The error bar shown at the physical 
point is only statistical and does not reflect any systematic uncertainties.) 



5.3 A first glance at the generalized isoscalar form factors of the nucleon 

In this section we present the results for the momentum and quark-mass dependence of the generalized 
form factors of the nucleon for the isoscalar flavour combination u + d at 0{p^) in BChPT. We note 
that at this order the only non-zero loop contributions to the three isoscalar form factors arise from 
diagrams c) and e) in Fig{Tl as the coupling of the isoscalar tensor field to the nucleon is not affected 
by chiral rotations. One obtains 



Cl3 



(51) 



with A2 0(0) given in Eq. (|44p and F|Q(t = 0) = 0. Interestingly, to the order we are working here, the 
t-dependence of this isoscalar form factor A2 Q{t) is given by the same function 



Floit) = Fl,{t) + 0{p' 



(52) 



that controls its isovector analogue Eq. (1411) , albeit with larger numerical prefactors (compare Eq. ()5ip 
to Eq. (j40|) ). We note that i^2o(0 does not depend on the scale A of dimensional regularization for 
the loop diagrams. The chiral coupling C13 is therefore also scale-independen10, parametrizing the 
(quark- mass independent) short-distance contributions to the radius/slope of ^2o(^)- "^^^ unknown 
contributions from higher orders in the chiral expansion can be estimated from a calculation of the 
triangle diagram displayed in FigJJJ Due to the coupling of the tensor field to the (long-range) pion 
cloud of the nucleon, among all the contributions at the next chiral order this diagram should give 
the most-important t-dependent correction to the covariant 0(j>'^) result of Eq. (j5ip . resulting in the 
estimate 0{p'^) ~ ^^ho (^'"^tt)- The explicit expression for this function can be found in appendix 



We note again, that this scale-independence refers to the UV-scales of the ChEFT calculation, not to be confused 
with the scale- and scheme-dependence of the quark-operators on the left hand side of Eq.((6]), which is completely outside 
the framework of ChEFT. 
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ID.2[ For completeness, we also note that in the limit 1 / {IGtt'^ Mq) ^ we obtain the corresponding 
HBChPT result 

^loitrHBChPT = alo + ^ml^ + ^t + 0{l/{l67r'F^Mo)), (53) 

which is just a string of tree level couplings. 

To order 0{p^) in BChPT the two other isoscalar form factors read 
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6mlM'^M'' + 6mlM^ ( - 3m^M^ + ) log 
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6ml 



4M2 



mt - 5mlM^ + hM'^ 



arccos 
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+ABUt,m^), 



(54) 



+6mlM^M'^ - 6mlM^ (m^ - 2M^'^ log ^ 

+ ^2S^= \mt - AmlM^ + 2mA arccos f 1 
^4M2 - ml ^ J V2M;/ 

+AC^.„Xt,m^). (55) 

with M defined in appendix O (cf. Eq. ()C.6|) ). MAr(m,r) denotes again the (quark-mass dependent) 
mass function of the nucleon Eq.(l9|), introduced via Eq.Q. In the limit l/(167r^F^Mo) — > we obtain 
the corresponding 0{p^) HBChPT results for i?|Q(t), C|Q(t), which at this order only consist of the 
tree- level couplings 61 O' '-2 o- As in the case of A2 ^{t) we have estimated the contributions from higher 
orders via 0{p^) ~ Ai?^ ^(t, m-,^), ^(^h o(^' "^tt), assuming that the dominant t-dependent higher order 
corrections to our covariant 0{p^) BChPT results of Eqs. (|54l55p originate from the 0{p^) triangle 
diagram displayed in FiglH Explicit expressions are given in appendix ID.2I We note that the non- 
analytic quark-mass dependent terms in 9(^)5 B2Q{t), C2Q(t) calculated in ref.[ll] with the help of 
the the HBChPT formalisnl^ correspond to the leading terms in a 1/{16tt^F^Mq) truncation of the 
0{p^) BChPT corrections AXh,o.{t, m^), X = A, B, C of appendix[El 

At this point we refrain from a detailed numerical analysis of the i-dependence of the generalized 
isoscalar form factors ^20(^5 -^2o(^)- hand very few lattice data for this flavour combi- 

nation have been published so far for pion masses below 600 MeV. Moreover, available lattice data 
neglect contributions from "disconnected diagrams" and are therefore accompanied by an unknown 
systematic uncertainty, which is very hard to estimate. On the other hand, in the t-dependence both 

^^In refs. [15l I16j additional terms have been calculated within the HBChPT approach. While some terms correspond 
to 0{p^) and C?(p^) contributions according to our power-counting, expanding the covariant 0{p^) result of Eq. (|54p to 

the order ^ leir^F^A/o ) '^'^^ ''^^Snize a term ~ 02,0"^^ present in ref.[15). However, as far as we can see, 

neither ref.[15) nor ref. |16) presents a complete 0{p'^) HBChPT calculation of the matrix element Eq.((7]). 
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of A|Q(t) and of -B|Q(t) we encounter chiral couplings (cis at OijP') in Eq. ([5T]) and -B34 at 0{p^) in 
^^ho (^j'^tt) of Eq. (|D.6p ) connected with (presently unknown) short-distance physics. The resulting 
leading i-dependence of these two generalized isoscalar form factors is hard to quantify at this point 
without any additional input either from experiment or from lattice QCD. We are therefore postpon- 
ing this discussion until further information is available, e.g. from the new lattice QCD simulations 
at low momentum transfer [39]. In the meantime we are preparing [32] a full (next-to-leading one 
loop) O^p^) BChPT analysis of the isoscalar moments of the GPDs which (in addition to several 
other diagrams!) also contains the contributions from the triangle diagram shown in Fig 121 already 
presented in appendix ID. 2[ 

Before finally proceeding to the summary of this work, we want to take a look at the third generalized 
isoscalar form factor C|Q(t) of Eq. ()55p . According to our power-counting, short distance contributions 
to the radius of this form factor are suppressed and only start to enter at 0{p^), both in HBChPT and 
in BChPT. After adding the 0{p^) estimate ACf^ ^ {t,m^) to the ©(p^) BChPT result of Eq.l^, we 
can hope to catch a first glance of the t-dependence of this elusive nucleon structure. Utilizing a| g of 
table 131 and assuming (x) 
its slope 



0.5 at a renormalization scale /x^ = 4 GeV^ [l3] we can determine 
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(56) 



At the physical point this would give us p^^mj^ = 0.14GeV) = —0.77 GeV ^. Truncating Eq. (|56p in 



l/(167r^F^A/o) we reproduce the chiral singularity 



~ m„ 



found in ref. 
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(57) 



It amounts to a slope of P(j\ 



-1.12 GeV 2, which is 45% larger than the BChPT estimate of 



Eq. (j56p . Interestingly, among the terms ~ ml shown in Eq. (j57p it is the 1 /Mo suppressed corrections 
to the leading HBChPT result of ref. |14j that dominate numerically. This gives a strong indication 
that a covariant calculation of AC^ ^ (t, m-j^) as given in appendix ID. 21 is advisable, automatically 
containing all associated (1/Mo)" corrections. 



^"Note that due to a different definition of tlie covariant derivative in tlie quark-operator on tlie left liand side of Eq.© 
our definition of tlie third isoscalar form factor differs from ref.[14| by a factor of 4: {t) = 4C2_o(i)- 
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The quark-mass dependence of the slope function of Eq. ([56|) akeady suggests that one obtains 
an interesting variation of the t-dependence of this form factor as a function of the quark-mass! We 
therefore close this discussion with a look at FiglHl There we have fixed the only unknown parameter 
C2 = —0.41 it 0.1 such that the BChPT result coincides with the dipole parametrization of the 
QCDSF collaboration at t = [22] for the lightest pion mass in the simulation, i.e. rriT^ = 640 
MeV. We note explicitly that this coupling only affects the overall normalization of this form factor, 
but does not impact its momentum-dependence. It is therefore quite remarkable to observe that the 
resulting t— dependence of this form factor according to this ChEFT estimate agrees quite well with 
the phenomenological dipole-parametrization of the QCDSF data at this large quark-mass, even over 
quite a long range in four-momentum transfer. The result is rather close to a straight lin^^^ (see Figj8]). 
We remind the reader that the value of this form factor at t = and rri-Tr = 0.14 GeV determines the 
strength of the so-called D-term of the nucleon, playing a decisive role in the analysis of Deeply Virtual 
Compton Scattering (DVCS) experiments [T]. Utilizing the extracted value of c| g we can now study 
the C-form factor also at the physical point, with the result also shown in FiglHl At this low value of 
the pion mass one can suddenly observe a non-linear t— dependence for low values of four-momentum 
transfer, due to the pion cloud of the nucleon. This is a very interesting observation, because such a 
mechanism would allow for a more negative value of C|Q(t = 0) at the physical point than previously 
extracted from lattice QCD analyses via dipole extrapolations {e.g. see ref. |22j). We obtain 



do {t = 0,m^ = 140 MeV) ^ -0.36 ± 0.1, (58) 

The assigned error corresponds to the fit error of C|Q(m^ = 0.64 GeV, t = 0)qcdsf given in ref. |22j . 
as it directly influences our unknown coupling c^q. However, we note that the (unknown) systematic 
uncertainties underlying the (quenched) simulation results of ref. [22] are not accounted for in this 
error bar. A comparison of these data with the new dynamical simulations of QCDSF and LHPC |39| 
might shed more light onto this unresolved question of systematic uncertainties in lattice 

With this value we can finally obtain the first estimate for the radius of this elusive form factor: 

f--^'' = Cj,,(t = 0,J.^140MeV) ^^''"- = "°"°^> 

^ (0.5±0.1)fm2. (59) 



We compare this result with the radii of the isovector Dirac and Pauli form factors of the nucleon, 
which are also dominated by pion cloud effects. Interestingly, with (r]')^ = 0.585 fm^ and (r|)^ = 0.80 
fm^ [U] the estimated value for (r^)^ seems to lie in the same order of magnitude! We note, however, 
that our numerical estimate of its slope Pq of Eq. (l56p given above is significantly smaller than the 
corresponding slopes of the isovector Pauli and Dirac form factors, as expected from general arguments 
and as already observed in lattice QCD simulations with dynamical fermions [45] (at very large quark 
masses) . 

However, before we can go into a more detailed numerical analysis of these interesting new form 
factors of the nucleon, one should first complete the 0{p^) calculation of the generalized isoscalar form 
factors, as there are additional diagrams next to Figl2l possibly also affecting the t-dependence, albeit 
presumably in a weaker fashion |32j . 

^^Most of the pion loop diagrams in BChPT with IR regularization as employed in appendix [Bl have the remarkable 
feature that they automatically go to zero in the limit of large pion masses, without any extra assumptions. This aspect 
of BChPT will be discussed in detail in ref. |18j. 

^^However, according to our understanding even the new dynamical simulation results of QCDSF and LHPC will 
neglect all contributions from disconnected diagrams. 
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Figure 8: A first glance at the momentum dependence of the form factor C|Q(t). The quasi-hnear 
t— dependence at the pion-mass m,r = 640 MeV (sohd hue) has been normahzed to the dipole 
parametrization of the QCDSF data of ref. [22] (dashed hue). The resulting non-linear t-dependence 
in this form factor for smaller values of is then due to the coupling of the tensor field to the pion 
cloud, providing an interesting mechanism to obtain "large" negative values at i = 0, 771,^ = 140 MeV. 




6 Summary 

The pertinent results of this analysis can be summarized as follows: 

1. We have constructed the effective chiral Lagrangean for symmetric, traceless tensor fields of 
positive parity up to O(p^) in the covariant framework of Baryon Chiral Perturbation Theory 
for 2 light quark flavours. 

2. Within this covariant framework we have calculated the generalized isovector and isoscalar form 
factors of the nucleon Ajol^'"^^)' -^2'o(*' "^^)' ^2o(^'"^^) to 0{p'^), which corresponds to 
leading-one- loop order. We can exactly reproduce the corresponding non-relativistic 0{p'^) re- 
sults previously obtained in Heavy Baryon ChPT by taking the limit l/(167r^i^^Mo) 0. Several 
HBChPT results published recently could not (yet) be reproduced, as they correspond to partial, 
non-relativistic results from the higher orders 0{p^ ,p^). 

3. According to our numerical analysis of the quark-mass dependence of the generalized form 
factors, we have noted that for i?2o'(^) C2Q{t) the observable quark-mass dependencies 
could be dominated by the (well-known) quark-mass dependence of the mass of the nucleon 
MNiml). This mass function appears in several places in the chiral results due to kinematical 
factors in the matrix-element used in the definition of the generalized form factors. Such a 
"trivial" but numerically significant effect is already known from the analysis of lattice QCD 
data for the Pauli form factors of the nucleon. 

4. The pion-cloud contributions to all three generalized isovector form factors at finite values of t are 
very small. The momentum dependence of these structures seems to be dominated by (presently) 
unknown short distance contributions. The situation in this isovector channel reminds us of an 
analogous role played by chiral dynamics in the isoscalar Dirac and Pauli form factors of the 
nucleon. At this point we were therefore not able to give predictions for the numerical size of the 
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radii or slopes of these interesting nucleon structure quantities. It is hoped that a global fit to 
new lattice QCD data at small pion masses and small values of t — extrapolated to the physical 
point with the help of the formulae presented in this work — will lead to first insights into this 
new field of baryon structure physics [39j . 

5. The BChPT leading-one-loop results for the generalized isoscalar form factors A2Q{t), -BfoCO 
are quite surprising. As far as the topology of possible Feynman diagrams is concerned, one 
is reminded of the isovector Dirac and Pauli form factors of the nucleon. A power-counting 
analysis, however, told us that those diagrams (e.g see Figj2]) from which one would expect 
large non-analytic contributions in the quark-mass only start to contribute at next-to-leading 
one-loop order. Our analysis therefore suggests that the momentum-dependence at low values 
of t is dominated by short-distance physics. A new quality of lattice QCD data for low pion 
masses and low values of t to constrain some of the unknown couplings is urgently needed before 
further conclusions can be drawn |39j . 

6. It is the value of {x),^ of a pion in the chiral limit that controls the magnitude of those long- 
distance pion-cloud effects in the generalized isoscalar form factors of the nucleon, pointing to 
the need of a simultaneous analysis of pion and nucleon structure on the lattice and in ChEFT. 

7. In the forward limit, the isovector form factor ^2o(^ ~^ ^) reduces to {x)u-d- Our covariant 
0{p^) BChPT result for this isovector moment provides a smooth chiral extrapolation function 
between the high values at large quark-masses from the LHPC collaboration and the lower value 
known from phenomenology. The required (chiral) curvature according to this new analysis does 
not originate from the chiral logarithm of the leading-non-analytic quark-mass dependence of 
this moment — as had been speculated in the literature for the past few years — but is due to 
an infinite tower of terms (m7r/Mo)" with well-constrained coefficients (see Eq. (j28p ). The well- 
known leading-one-loop HBChPT result for (a;)„„rf of Eq. (j30p was found to be applicable only 
for chiral extrapolations from the chiral limit to values slightly above the physical pion mass, as 
expected. 

8. Judging from the available (quenched) lattice data of the QCDSF collaboration our O(p^) 
BChPT result of Eq. (fH]) for ^|o(^ = 0) = {^)u+d also provides a very stable chiral extrap- 
olation function out to quite large values of effective pion masses. 

9. A study of the forward limit in the isoscalar sector has led to a first estimate of the contribution 
of the u— and d— quarks to the total spin of a nucleon Ju+d ~ 0.24. This low value compared to 
previous determinations arises from the possibility of a small negative contribution of -B|Q(t = 
0) ~ —0.06 at the physical point, driven by pion cloud effects. However, at the moment the 
uncertainty in such a determination is rather large, due to the poor situation of available data 
from lattice QCD. 

10. In a first glance at the third generalized isoscalar form factor C|g(t) the quark-mass dependence 
was found to be qualitatively different from A2 o{t), SI o(^)- slope contains a chiral singularity 
~ m^j^^ and the influence of short distance contributions is suppressed. A first numerical estimate 
of its slope gives f« —0.75 GeV^, which is much smaller than the slopes of corresponding 
Dirac or Pauli form factors. At low t we have also observed signiflcant changes in the momentum 
dependence of this form factor as a function of the quark-mass, resulting in the estimate C"! o(i = 
0) ~ —0.35 at the physical point. 

11. Throughout this work we have indicated how to estimate possible corrections of higher orders 
to our leading-one-loop BChPT results. The associated theoretical uncertainties of our 0{p'^) 
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calculation have been discussed in detail. Ultimately, in order to judge the stability of our 
results it is mandatory that we analyse the complete next-to-leading one-loop order. Work in 
this direction has already started [32]. 

Finally, we note that the tensor Lagrangeans constructed in section [3] invite a host of further studies, 
pertaining both to generalized axial form factors of the nucleon [29] and to the energy-momentum- 
tensor of the nucleon [33] ■ 
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A Basic Integrals 



The integrals required for one-loop calculations in BChPT can be reduced to two basis integrals in 
d-dimensions: 
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where m (M) is a mass function involving the mass of the Goldstone Boson (of the Baryon) and p^^ 
denotes a four-momentum determined by the kinematics. The propagators are shifted into the complex 
energy-plane by a small amount e to ensure causality. Utilizing the MS-renormalization scheme of 
ref.[l] with 
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one obtains the dimensionally-regularized results [1] 
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More complicated integral expressions needed during the calculation are defined via 
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The integrals H^l are related to the two basis integrals of Eqs. ()A.lllA.2p via well-known tensor- 
identities [U [18] . Finally, we note that the integrals involving more than one baryon propagator can 
be related to the ones defined above via 



(A.8) 
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B Regulator Functions 

The basis regulator function is defined as \17\ [T8] 

Rii[M\m\p') = j dxj j^[xM^ + {a 



x)p^ + (1 — x)m^ — / 
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More complicated regulator functions R 
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11 ' 



1 ... 5 can be defined in analogy to Eqs. ()A.6|A.7p . We 



note that for our 0{p'^) calculation of the moments of the GPDs we only need to know these functions 
up to the powej^ of m^, t in order to obtain a properly renormalized, scale-independent result, which 
at the same time is also consistent with the requirements of power-counting. The regulator functions 
needed for our 0{p^) BChPT calculation (see Appendix [C]) read 
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The derivatives of the regulator functions needed for the calculation (see Appendix [C]) read 
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One can clearly observe that all contributions are polynomial in m2 (and therefore polynomial in the 
quark-mass jl7j ) or polynomial in t [18], as expected. Their addition to the MS-results therefore just 
amounts to a shift in the coupling constants |17[ [15] of the effective field theory and does not affect 
the non-analytic quark-mass dependencies, which are the scheme-independent signatures of chiral 
dynamics. 



■^■'strictly speaking we need to know the regulator terms contributing to the generalized form factor ^2'o(*) up to 
power {m^,ty, whereas for _B2'o(t), C2o{t) only the leading terms {m^,t)'^ are required [18] . 
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C Isovector Amplitudes in 0{p'^) BChPT 

The five 0{p'^) amplitudes in the isovector channel corresponding to the five diagrams of FiglUwritten 
in terms of the basic integrals 



of Appendices El and |B] read 
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Note that the various couplings and parameters are defined in section [31 r] denotes the isospin doublet 
of proton and neutron. The variables in the integral functions are given as 
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where t = corresponds to the momentum transfer by the tensor fields. 

Zj^f denotes the Z-factor of the nucleon, calculated to the required 0{p'^) accuracy in BChPT. It is 
obtained from the self-energy Sjv at this order via the prescription 



2^ = 1 + 



as 



N 



f=Mo 



+ 0{p% 



with 



-'N 



4F| 



(Mo +p)[mlh^{Mlml,p') + (Mq - (M^, _ a^] + 0{p 



(C.7) 



(C. 
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D BChPT Results in the Isoscalar Channel 



D.l Isoscalar Amplitudes in Oip"^) BChPT 

To 0{p'^) in BChPT the results in the isoscalar channel are quite simple. The amplitudes corresponding 
to the Feynman diagrams of Fig(T] can be be simply expressed in terms of results already obtained in 
the isovector channel discussed in the previous section [Cj They read 



Amp^, 



Amp^ 

Ampd 
Amp^ 



+ 



— Amp" + 0{p'^ 



Amp'' + 0{p-^). 



(D.l) 
(D.2) 
(D.3) 
(D.4) 



Note that the various couplings and parameters are defined in section [3l 



D.2 Estimate of 0{p^) contributions 

The contributions from the (higher order) Feynman diagram shown in figure [2] to the generalized 
isoscalar form factors read 



AAl^Xt,m^) 



2 du 



,9 2 p / -8 1 (3"^' + ^0 ) - + (2M2 + m^) / 
6Z7r j_i p I 

2M^ {IIM^ - Im^M^ + 2m^) p'^ + 12Mo^ [M^ - m^fp^ 

/ + 2,f (3Mo^ - m'^) - 2p^ {M^ - m^) {7M^ - Irh^M^ + m^) 

+6Mo2(Mo2-m2)^]log^ + 



2M2 



2M2 (p2 + ^2) _ (p2 _ ^2)2 _ 

6Mo^° - 4M^ (Gm^ + 5p^) + (36m^ + 38mV + 23/) 



-M^ (24m^ + ISp^m'^ + llp^m^ + 9/) + {6m^ + 2p'^m^ 



-bp'^m^ - 2fm^ - f) + f (f - m^ f (p^ + 2m^) 



arccos 



Ml + rr? 



2Mom 



(D.5) 
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967r2F2 Mn 



du 



Mi 



6A/n + 18m: 21og — ^ - 3 + t 11 - 61og — ^ 



A 



Mn 



A 



36Mo^ 



- (4m2 + 3p2) 
2/ - 3 {2,M^ - m^) 



Q fi '2 

+771 — m p 



arccos 



2Mom y p' 

6 (Mq^ - m2)^p2 ^ 6 (MqV - 2Mo^ (Mq^ - ?^^) + (Mq^ - ?^^)^) log ^ 



33 V / 71 //^ A 2 34 V / 



Mo A2 



Mo A2' 
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+9Mo^ 



(^-mV + 2m V {rh'^ - M^) + (Mo^ - m^f^ 
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log 
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Mo 



9M2 
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-4 (4m2 + p2) ^^2^6 _ (_24n2j^4 _^ 4y2^2j^2 _^ ^4^ ^4 

+ {-IGu^m^ + 20uV?Ti^ + 2 (1 - 2u^) ih^f + /) Mo^ 



+m^ (m^ - ijku^frv' - Su^m^p^ + {iu^ - l) p'^) 



arccos 



M^ + m2 - ^2 
2Mom 



with the new variable 



m~ = rrii + { u'^ — - \ t. 



(D.7) 



(D.8) 



We note that the contributions of AA2o(^''^t)' ^C*! 0(^5 ^^17^) are finite at ©(p^) in BChPT, while 
AB2 Q{t,mT^) contains two new counterterms B^^^{X), 534(A) at this order |32j . 
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